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1. AJhatracr 


Future robotic manipulators carried by a spacecraft will be required to perform complex 
tasks in space, like repairing satellites. Such applications of robotic manipulators will encounter 
a number of kinematic, dynamic and control problems due to the dynamic coupling between 
the manipulators and the spacecraft. This p a per p re s e n te a, new ana lytical modeling method 
for studying the kinematics and dynamics of manipulators in space. /'The problem is treated by 
introducing the concept of a Virtual Manipulator ( VM). The kinematic and dynamic motions of 
the manipulator, vehicle and payload, can be described relatively easily in terms of the Virtual 
Manipulator movements, which have a fixed base in inertial space at a point called a Virtual 
Ground. It is anticipated that the approach described in this pipsi-will aid in the design and 
development of future space manipulator systems. 


- / 


2, Introduction 


Robotic manipulators are potentially very useful for performing complex tasks in non-industrial hostile environ- 
ments [1,2|, such as in space. A number of studies have considered the potential applications of manipulators is 
space and the capabilities that these systems must have to achieve anticipated mission goals [3-5). These applications 
include tasks such as repairing, servicing and constructing space stations in orbit. Currently, these tasks can only 
be performed by astronaut Extra Vehicular Activity (EVA). Eliminating the need for EVA would obviously reduce 
hazards to the astronauts and mission costs. 


Unfortunately, the use of manipulators in space is complicated by the man ipulator/spacec raft dynamic coupling. 
For example, movements of a manipulator will disturb the attitude of the spacecraft carrying it. This coupling 
will adversely affect the manipulator’s precision, and reduce the on orbit life of the system by consuming excessive 
attitude control fuel. Also, any motions of the spacecraft, say due to the firing of attitute control jets, will disturb 
the manipulator. Therefore, new manipulator concepts, designs and control techniques will be required to minimize 
and compensate for the manipulator/spacecraft dynamic coupling. 

Researchers working on the control of space manipulators have focused their attention on issues such as sen- 
sor reqirements, path planning algorithms, teleoperator control (6-8] ; the problems of vehicle/manipulator dynamic 
interactions remain unresolved. 


This paper presents a new and effective analytical modeling method for studying the kinematics and dynamics of 
manipulators in space. The problem is treated by introducing the concepts of a Virtual Ground (VG) and Virtual 
Manipulator (VM). As discussed below the VG is located at the center of mass of the manipulator/spacecraft system, 
and the VM is an ideal kinematic chain connecting the VG to any point on the real manipulator. Motions of & 
system, including a vehicle, manipulator and payload can be described easily by the VM. This model has proven to 
be effective in calculating the kinematic and dynamic properties of the system; such as its inverse kinematic solution 
and workspaces. This paper shows that the VM approach can also be used*to plan the manipulator’s motions in 
order to minimize the degrading consequences of the manipulator/spacecraft dynamic interactions. 


3, A Model of Manipulators in Space 

Future space manipulator systems will have one or more mechanical arms carried by a vehicle, as shown in 
Figure 1. The vehicle will be capable of motion in six degrees of freedom, and will have reaction jets for position and 
attitude control. Although manipulators could be driven by photovoltaicly powered electric actuators, which use no 
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where 


( 2 ) 


N 

Mm = £ 

7 3 1 

Since there are no external forces, the VG is stationary in the frame N and the vector V § is always constant. 

In the following development the VM properties such as link dimensions and joint axes, for initial manipulator 
configuration axe described. Then the rules for its joint movements as a function of the real manipulator joint 
movements are presented. Referring to Figure 3, which shows the end effector VM for the manipulator shown in 
Figure 2, the i tk link of the Virtual Manipulator is defined by the vector V,*, 

Vx = D! 

V, = Hi+D, 

V N = H*-,+D* (3) 

where 

Di = Ri £ M,/M w (.• =1,2 N) (4) 

and f»i 

H, = U £ M' /Mm ( * =1,2 ft- 1) (6) 

f*l 

The first VM link represents the vehicle’s orientation. This link is attached to the VG by a spherical joint and 
its motion is equal to the three vehicle rotations with respect to inertial space. The end of the Virtual Manipulator 
terminates at the end effector, defined by a vector E, fixed in the N th VM link. 

The i th VM joint is taken as a revolute or a prismatic joint depending upon whether the i tk joint of the real 
manipulator is revolute or prismatic. The axis of rotation for a revolute VM joint, j4, is parallel to the axis of the real 
manipulator joint Similarly, the translational axes of prismatic VM joints are parallel to the corresponding axes 
of the real manipulator prismatic joints. Equations (1) through (5) define the VM and its position corresponding to 
the initial position of the system, as shown in Figure 2. The VM links will all be parallel to the real manipulator 
links in cases where all the centers of mass for all manipulator links lie on a line connecting the manipulator joints 
on the corresponding link. 

The VM will move as the joints of the real manipulator move. The angular rotations of the VM revolute joints, from 
their initial position, are equal to the angular rotations of the corresponding revolute joints for the real manipulator. 
The prismatic virtual joint translations are ratios of the corresponding real prismatic joint translations. For an end 
effector VM, translation of the virtual joint, P,, is given by: 

P i =T i t,MjM M ( 6 ) 

*=i 

For the VM in its position of construction, its initial position, the values of are taken as zero. Hence the initial 
magnitudes of Pj are zero. The prismatic joint motions, T, , are referenced to the initial position. 

If a VM that is constructed according to Equations (1) through (5), moves with the real manipulator according 
to the above description, and its link shapes and lengths remain constant as the manipulator moves, then it can be 
shown (see Appendix A) that: 

1 The axis of the i tK virtual joint is always parallel to the i th axis of the real system joint, and 

2 The Virtual Manipulator end point will always coincide with the real manipulator’s end effector. 


These properties enable the kinematic and dynamic motions of a free-floating manipulator system to be described 
by the motions of a much simpler Virtual Manipulator which has a fixed base in inertial space. The properties of the 
VM remain the same as long as the mass property of the system does not change. For example, when the manipulator 
grasps a free-floating pay load, the VM changes. According to Equations (1) through (5), the VM link lengths will 
be reduced for the addition of a payload. Virtual Manipulators constructed for points other than the end effector 
have different links than the links defined in Equation (1); and their joint movements maybe different than the ones 
described above, for example, prismatic joint translations may be the vector (P< - T<), depending upon location of 
the point used to construct the VM [9|. 
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reaction fuel, manipulator motions could disturb a vehicle's position and attitude and result in the consumption of 
excessive amounts of attitude fuel. The useful ! ! fe of spacecraft systems is often limited by the amount of reaction 
jet fuel they can carry. * 

Two approaches to solve this problem are: 1) permit the vehicle to move and compensate for the base motions in 
the manipulator task planning; and 2) plan the manipulator motions so that they do not cause the vehicle to move 
excessively. The first of these approaches requires the ability to perform inverse kinematic and workspace calculations 
for a free-floating system [10|. The second requires methods for planning manipulator motions that would self-correct 
the vehicle's orientation with little or no reaction jet adjustments. These approaches and associated issues are 
addressed here through the Virtual Manipulator technique. Assumed in this work is that the external forces/torques 
acting on the system are negligible, and that the system is free floating. Also assumed is that the system elements 
may be modeled as rigid bodies. The later assumption may not be valid if a manipulator must perform high speed 
motions. 

4* Analytical Development of the Virtual Manipulator 

The Virtual Manipulator (VM) is a massless kinematic chain terminating at an arbitrary point on the real 
manipulator. Its base is the Virtual Ground (VG), which is an imaginary fixed point in inertial space. It is proven 
below that for a given system the properties of the VM and location of the VG are fixed. VMs exist for many 
different manipulator structures, such as open or closed chains, single or many branch arms, revolute or prismatic 
joints [9-1 ij. The discussion in this paper will be limited to manipulators composed of spatial serial chains with 
revolute or prismatic connections. Although VMs exist for any point on the real manipulator, this paper deals with 
VMs whose end points coincide with the real manipulator end efTector. 

The VG is defined to be the center of mass of the manipulator system. From elementry mechanics, when there 
are no external forces, such as from reaction jets, the VG will be fixed in an inertial space. It will not move due to 
any internal forces of the system such as joint torques, or due to any manipulator motions. 

Figure 2 shows a schematic drawing of an N body spatial manipulator system. The first body in the chain 
represents the vehicle which carries the manipulator. The N th body is a combination of the payload and the last 
link. The joint is called Ji , and is the center of mass of the i th body. The vectors R* and L, connect C% to 
Ji and Ji to C, + i, respectively. The vector Ryv connects C * to the end effector. The vectors R* and Lj_| are fixed 
relative to the i th link, and hence the angle between these vectors is constant for all system configurations. If the t 1 * 
manipulator joint is a revolute joint, the vector defining the axis of rotation of J, is called A,, and the angle £, is the 
rotation of the t th joint. If the i th manipulator joint is a prismatic joint, the vector T < is defined to be the translation 
along the translational axis. If the i tH joint is revolute, then the magnitude T, is equal to zero. 




The location of the VG for this system in inertial space, the center of mass of the system, can be found by knowing 
some initial position of the system. The vector S(0) defines the initial known location of the end effector with respect 
to an inertial reference frame N. Then the location of the VG, the vector V tf , can be obtained from conservation of 
linear momentum by: 

N i-l 

V, = D s (°) - B R > + L > + T,)| MJM," (1) 

.=1 ;=l 
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Figure 3: N Body System and iU VM. 



Figure 4: A Three Body Planar System and 
iU VM. 


Table 1 gives the properties of a very simple planar manipulator and its Virtual Manipulator, shown in Figure 4. 
It should be remembered that the method is not restricted to planar systems. 


5. Applications of Virtual Manipulators 

The Virtual Manipulator approach has a number of possible applications. VMs can be used to simplify the 
inverse kinematics of space manipulators, calculate their workspaces, plan their motions and formulate the equations 
of motion [9-11]. It should be noted that ut’ng conventional methods, these problems are far more difficult for space 
manipulators than for industrial manipulators with fixed bases. In the sections below, the use of the VM is shown 
for workspace analysis and path planning. 

A. Workspace Analysis 

Since the vehicle and manipulator'dynamics are coupled, the manipulator’s motions will cause the vehicle to move 
and this in turn makes it difficult to find the manipulator workspace. In fact several different types of workspaces 
exist. In this section, a workspace called the constrained workspace, for a manipulator in space is defined, for a more 
complete discussion of space manipulator workspaces refer to references [9,10], For the constrained workspace it is 
assumed that the attitude, but not the location, of the vehicle is controlled. This can be achieved without the use of 
attitude control fuel by employing reaction wheels, or by using the self correcting maneuvers discussed later in this 
paper. 

To find the constrained workspace, a Virtual Manipulator is constructed to the end effector of the real manipu- 
lator. The joint limits of the real manipulator are transformed into VM joint limits. The workspace of the Virtual 
Manipulator is then found using conventional workspace analysis methods (12|. The real manipulator workspace will 
be equal to the VM workspace because of the following reasons. The VM end point coincides with the real manip- 
ulator end effector, and it is assumed that it is possible to control the orientation of the Brst VM link, representing 
the vehicle, with respect to inertial space. The other joints are controlled with their actuators. This workspace will 
always be a spherical shell, assuming there are no limits on the vehicle orientations. Figure 5 shows the constrained 
workspace for the simple two link manipulator shown in Figure 4, it was found using its VM. 

B. Path Planning 

In certain cases, the magnitude of the rotations of the vehicle caused by the manipulator’s motion may not be 
acceptable. For example, vehicle rotations may cause communication devices to loose their signals. Vehicle rotations 
can be controlled using reaction wheels or reaction jets. However, these devices have the disadvantages of increased 
mechanical complexity and system weight or increased consumption of attitude control fuel. 

It is shown below that the manipulator itself can be moved in such a way as to have the end effector follow a 
nominal specified path and yet have a prescribed vehicle orientation, with specified limits, without using attitude 
control fuel or requiring reaction wheels. 
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Tkble 1: Characteristics of Planar Manipulator 
and ita VM. 
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Figure 5: Constrained Workspace of Manipulator 
Shown in Figure 4. 

To find this motion the principle of conservation of angular momentum is applied to the system. For an n degrees 
of freedom space manipulator the following equation can be writter. 

X = P(e,X)A, (7) 

F: 3 by n matrix with elements 

X: 3 by 1 vector of vehicle inertial orientations with elements, X,, 

0 : n by 1 vector of joint angles with elements #, 


In general, Equation (7) is non-integratable, that is: 


d0 k d9; 


d*x t 

dijdiu 


Therefore, the final vehicle orientation depends on path taken by the manipulator from one position to another. It 
follows that the final vehicle orientation will change if the manipulator moves along one path in joint space sad 
returns to its initial position by another path. This is a similar notion to the one which permits astronauts to reorient 
their bodies by moving their limbs (13). This leads to a strategy for adjusting or correcting motions of the vehicle's 
orientation. In this strategy nominal trajectories are selected for the end effector and vehicle orientation. Then the 
joint motions are executed assuming the vehicle follows its trajectory. If at any point the vehicle orientation deviates 
from its desired path by more than a specific amount, a series of small cyclic motions, selected to correct for the 
vehicle orientation are added to the joint motions. 

To find the cyclical joint motions that achieve the desired vehicle/ baae orientation corrections, it is assumed that 
these motions are small enough that the end effector deviates only by a small amount from it# nominal trajectory. 
This small motion assumption permits the use of a nonlinear system model in which nonlinearities of order greater 
than 2 can be neglected. 

First, let X be a set of Euler angles defining the baae orientation with respect to an inertial coordinate frame. 
The initial and desired final baae orientations are X« and X#, respectively. The desired change in the Euler angles is 
defined by 

SX = X, - X d (9) 

Let 0o be the vector defining the initial and final joint positions at the beginning and end of the correction 
maneuver. Also let the vectors 6V and £W define small joint movements. The cloeed correction path is constructed 
by having the manipulator move along the straight lines, in joint space defined by vectors £V and 6W shown in 
Figure 6. 

For small 5V and 6 W the following equation can be obtained from Equation (7). 


3 3 3 

1 ; = t m= 1 




(*= 1.2.3) 


( 10 ) 




where 6X it 5VJ and BWi are elements of the vectors 6X t 6V and $W, respectively. In the case of a three DOF spatial 
manipulator, Equation (10) will yield three equations with six unknowns. Three additional constrain equations are 
required to solve for SV and S W. 

If vectors 6V and 6 W are parallel, the cyclic motion will not produce any vehicle rotation. Therefore it is assumed 
that these vectors are perpendicular: 

6V t . *W = 0. (II) 

Further, the magnitudes of 6W and 6V are assumed to be equal: 

6\ t • 6\ = 6 W r • $W, (13) 

and one of the elements of is chosen to be a linear combination of the other two. For example, 

*V, = (*V, + *Vi)/2 (13) 

Equations (10) through (13) yield six scalar equations with six scalar unknowns, which can be solved for the 
desired joint trsjectories, 6V and 6W. If the required correction, 6X, is large, the values of 6V and 6W may violate 
the small joint motion assumption. In this case the desired correction can be achieved by a series of m cyclical 
correction maneuvers. It is shown below that at each cycle Equations (10) through (13) do not have to be resolved 
and the final position can be achieved. 

Referring to Figure 7, T(X ; ) is a 3 by 3 matrix which transforms a vector expressed in vehicle body coordinates 
(x,y,i) into inertial or Newtonian coordinates (S 9t S 9t N M ) t when the body is at jth orientation. The transformation 
matrix for the initial vehicle orientation is T(X*). The transformation matrix for tbs desired vehicle position to be 
achieved after m cycles is T(X^). After one correction cycle, the transformation matrix ia T(X* + 6x) t where, 

T(X. + *x) = T(X,)A, (14) 

and the matrix A is the transformation matrix from the vehicle position, one cycle from the initial vehicle position, 
back to the initial position. The A matrix will not change with each cycle because the total system, vehicle asd 
manipulator, have been subject only to a rigid body rotation in inertial space. Hence after m cycles the transformation 
matrix from the desired system position to inertial coordinates is simply: 

T(X^) = T(Xt)A m (15) 

Equation (15) can be solved for A: 

A = PA , /’ n p-' (1«) 

where A is a diagonal matrix of the eigenvalues of T(X, ) ” 'T(X^) and P is a matrix of corresponding eigenvectors. 

Using the A matrix obtained from Equation (16), the change in Euler angles ($x) are calculated from Equa- 
tion (14). Then the joint correction motions for each cycle, <5V and £W, are obtained by solving Equations (10) 
through (13) The manipulator should go through the derived joint transformations (6V, <5 W ) m times to approach 
the desired vehicle orientation. However, the final vehicle orientation after m cycles will usually be slightly different 
than the desired orientation because of the neglected higher order nonlinearities. In order to achieve the desired 
vehicle orientation more precisely, the over all correction may need to be broken into several smaller corrections and 
the process repeated with a slightly different set of 6V and 8 W for each subcorrection. 




Figure 6: A Closed Path Correction in Joint Space Figure 7: Vehicle Coordinate Rotation Due to Cyclic 
for a Vehicle Rotation. Manipulator Motion. 
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The above technique is now demonstrated for a spatial 3 DOF space manipulator shown in Figure 8. The 
properties of the manipulator are given in Table 2. It is desired to rotate the vehicle from its initial orientation to 
its Snal orientation as shewn in Thble 3. hi this example, it was necessary to solve for the joint trajectories (3V and 
IW) 3 times to precisely obtain the desired vehicle orientation. The joint tr^ectoriss for these 3 cycle sets are shown 
in Figures 0 through II. Bach cycle is repeated 30 times to achieve the desired vehicle orientation. Thble 3 shows the 
system angles after each cycle set. During each cycle the vehicle os cill a t es in sympathy to the manipulator's motion, 
see Figure 12. However the mean orientation of the vehicle changes continuously and reaches at the same time 
the joints return to their initial positions. Figure 13 shows the mean vehicle orientation during the joint cycles. The 
vehicle movements during the joint motions are ± 0.1 radians from their mean trajectory. Here one can clearly see 
change in the base orientation as the manipulator joints cycle through their motion. 
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Figure S: Spatial 3 DOF Space Manipulator. 
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Figure 10: Joint Angle Trajectories for Second 
Set of Cycles. 



Figure 9: Joint Angle Trajectories for First Set 


of Cycles. 



lhble 3: Manipulator Angles at Different Instances. 
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Table 2: Three DOF Manipulator Characteristics. 
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Figure 12: The Xj Vehicle Coordin.te for the Fi « UM 13: Me «» Vehiel « Eulw During 

First Set of Cycles. Correction Cycle. 

0. Summary and Conclusion 

In this paper, the concepts of Virtual Manipulators and Virtual Grounds are discussed. The end effector VM 
characteristics and proof of its properties for serial link with revolute and prismatic joints were presented, and some 
of its applications were discussed. This is a new concept and further research is required to demonstrate its full 
capabilities. 
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Appendix At Proof of Virtual Manipulator Properties 

First it will be proven that for a VM constructed using the rules presented in section 4 of this paper the VM end 
point will coincide with the end effector. Then it will be proven that when the manipulator goes through a movement 
the VM joint motions described in section 4 will keep the VM end point on the end effector. 

First, recognising that the system center of mass is stationary in the inertial frame N, V # remains stationary in 
this frame and referring to Figure 2 yields: 


N-l 

M, 0 ,V, = M n 8 + - L s-i - Tw-i - R*-,| + . . . + M,(S - £ (U + T< + R,)| (Al) 

Recall that if the joint is revolute T< = 0, otherwise, its magnitude is equal to the prismatic joint translations 
from the initial manipulator configuration and its direction is along the translational axis. Equation (Al) can be 
solved for S(t) as follows, 

S W . V, + + *, + 1.) + . . . + «■ + »»♦■ +«»-' ( R„_, +Tj,_i + Lff-i) (AS) 

Mtoi Mtot 

Equation (A2) can be written in terms of the vectors D # , H, and P. by using Equations (4) through (6), to yield: 

S(t) = V, + (Di + Hi + P t ) + . . . + (Dur-t + H*-i + Pjv-i) (A3) 

Using Equation set (3) and the fact that the end effector position is always equal to S(t) + R/v gives: 


E(f) = V, + V! + Pj + . . . 4- Ptf-i + V* (A4) 

It should be noted that this equation does not depend upon the existance of the Virtual Manipulator. The vector 
chain represented by Equation (A4) describes the end effector position relative to the N reference frame for all time. 

For the initial manipulator position the VM constucted according to the procedure outlined in section 4 has an 
end point described by the following vector chain, 

V, + V 1 + ... + V* (A5) 

Comparing Equations (A4) and (A5) it follows that in the initial position, when the P/a are equal to zero, the end 
effector coincides with the VM end point. 

Now it will be proven that as the real manipulator moves the VM joint motions described in section 4 will keep 
the VM end point on the real end effector. Say the manipulator goes through some joint movement, from section 4, 
the following vector chain describes the VM end point, where the P,’s are no longer zero, 

v; + v; + p; + . . . + p*_ t + (A 6) 

In the following paragraphs it will be proven that the vectors V?, V * and P* in Equation (A6) are the same as V*, 
V, and P f in Equation (A4), respectively, therefore, the VM end point will coincide with the real end effector. The 
vector V # is always constant, therefore, 

v ;=v, 
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(A7) 


The initial real manipulator links are composed of vectors Lj.j + R, and since the manipulator links are rigid, the 
magnitude of the vectors L,-j and R< and the angles between them are always constant. Since the magnitudes of 
and R< are constant, then from Equations (4) and (5) the magnitudes of and D< will also be constants. It 
can also be seen that the angles between H«-i and D, are constant. Then 

|V,*| = |H*-i + D«| VI. ■ (A8) 

The Virtual Manipulator links are composed of the Vf vectors. These links don't change their shapes and lengths 
as a function of time and since magnitudes of Vf are initially equal to magnitudes of V*, and magnitudes of V, do 
not change with time it follows that: 

|V;| = |H,_ 1 + D < | = |V < | Vt,.’ (A9) 

The magnitude of the vectors Pj in Equation (A4) and the PJ vectors in Equation (A6) are both obtained from 
the real manipulator prismatic joint translations, using Equation (6), therefore by definition, 

|p;i = |P<l *,.• (aio) 

Now it will be proven that the direction of the vectors Vf and P t * in Equation (A6) are parallel to vectors V,* and 
P< in Equation (A4), respectively. 

First it can be established that the rotations of the first VM link are set equal to the vehicle rotations and hence 
the first VM link will always be parallel to the vehicle. Therefore, 

V;=V! Vt (All) 

Since axis of rotation or translation of the first real and virtual joints are fixed relative to their corresponding first 
links, and the rotations of the first VM link is the same as the vehicle, and these axes are initially constructed to be 
parallel, then they will always be parallel. 

Now consider the case when the first real manipulator joint is revolute. The elements of the second VM link H{ 
and DJ will be parallel to L| and Rj and in turn the vectors VJ and Vj will be parallel and 

VJ = V, Vt (A12) 

because the rotational axis of the first VM and real manipulator joints are parallel, as shown above, and the magnitude 

of their rotations are equal by construction. 

In the cases where the first joint is prismatic, the elements of the second VM link HJ and D$ will be parallel to 
L| and R* and in turn the vectors VJ and V 3 will be parallel and 


VJ = V* Vt (A13) 

because the VM and real manipulator translational axis for this joint are parallel. In the same manner it is possible 
to show that 

v; = v< vt,,- (a 14) 

Also, in a similar manner all the translational axis of the real manipulator and the VM are always parallel, and from 
Equation (All), 

P*=P< Vt,,- (A15) 

Substituting Equations (A15), (A14) and (A7) into (A6), and comparing the result with Equation (A4) shows that 
the end effector will always coincide with the VM end point, and this completes the proof. 
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